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Quasiclassical Theory for Thermodynamic Properties of the Two-Band 
Superconductor MgBg under Magnetic Fields 

Masatoshi MuGiKURA*, Hiroaki Kusunose and Yoshio Kuramoto 

Department of Physics, Tohoku University, Sendai 980-8578 

Consistent account is provided of strongly anisotropic properties in thermodynamics of MgBj 
under magnetic field H using approximate analytic solution in the quasiclassical formalism. 
With the strength of the pairing interactions and the Fermi velocities of two bands as fitting 
parameters, it is possible to reproduce in wide regions of the phase diagram the H- and tem- 
perature dependences of relevant quantities such as the zero-energy density of states (ZEDOS), 
two energy gaps, the specific heat and the Hc2 anisotropy. It is demonstrated that the velocity 
ratio of two bands affects considerably the low-field behavior of the ZEDOS and energy gaps. 
The ratio also affects the high-temperature behavior of the Hc2 anisotropy. The anisotropy of 
the coherence length in the a band depends weakly on H and T, whereas the coherence length 
in the vr band remains isotropic even in the presence of the inter-band coupling. 
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1. Introduction 

Superconductivity in MgBj^ has attracted much inter- 
est not only for its high transition temperature Tc = 39K 
but for the unique character coming from distinct Fermi 
surfaces in a and tt bands. Owing to the presence of 
multiple bands, the superconducting state has two differ- 
ent energy gaps, although the pairing is triggered by the 
electron-phonon coupling in the a band. The two-band 
nature of superconductivity has been supported by vari- 
ous experiments such as point contact spectroscopy,^''^ 
specific heat,''"^ and the angle-resolved photoemission 
spectroscopy (ARPES).^ In addition to providing two 
different gaps, the a and tt bands have different dimen- 
sionality, i.e., the former is quasi two-dimensional and 
the latter is rather three-dimensional. The interplay of 
the two gaps with different size and anisotropy should 
lead to characterisitic behaviors, particularly in thermo- 
dynamic quantities in a mixed state. Indeed, the strong 
temperature dependence of the anisotropy of the upper 



critical field, r(r) = H^^f\T)/H^^^'{T), has been ob- 
served*"^*^ and has been discussed theoretically as a con- 
sequence of such interplay.^^"^^ 

There have been a number of theoretical investiga- 
tions on various thermodynamic properties in the mixed 
state, such as the zero-temperature specific-heat coeffi- 
cient 7(_ff)^^'^^ and the field dependence of the energy 
gaps.^^"^° However, validity of each theoretical tool used 
in these studies is limited to either H ~ Hc2 or T ^ 0. 
Hence consistent account of the whole H-T phase dia- 
gram within the single framework is highly desired. 

Recently, one of the present authors has developed 
a theoretical tool to investigate thermodynamic prop- 
erties^^ using the approximate analytic solution in the 
quasiclassical formalism.^'* The approximate solution 
was originally proposed by Brandt, Pesch and Tewordt 
(BPT) in the Gor'kov formalism, and is most suitable 
to strongly type-II superconductors near Hc2- However, 
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a comparison with reliable numerical calculations indi- 
cates that the present approximation is applicable semi- 
quantitatively to the wide region of the H-T phase di- 
agram. The purpose of this paper is to give a consis- 
tent account of thermodynamic properties based on the 
approximate analytic solution, especially of anisotropic 
behaviors in the mixed phase. 

The paper is organized as follows. In the next section, 
we briefly introduce the method of approximate analytic 
solution. In § 3, we first determine the plausible strength 
of pairing interactions at = 0, and discuss the ratio 
of the Fermi velocities of the two bands focusing on the 
temperature dependence of the anisotropy of Hc2 ■ Then, 
using the parameters obtained, we compare H and T de- 
pendence of the specific heat with experimental results. 
It is shown that the \ow-H behavior of the zero-energy 
density of states (ZEDOS) and the gaps is sensitive to 
the ratio of the velocities of two bands. The H and T 
dependence of the anisotropy of the coherence length is 
also discussed. The last section summarizes the paper. 

2. The BPT approximation 

To discuss thermodynamics under magnetic fields we 
have to solve the Eilenberger equation, the Maxwell 
equation and the gap equation self-consistently. It is te- 
dious to solve these equations as they stand. However, 
if we approximate the spatial dependence of the internal 
field and the diagonal part of the quasiclassical Green 
function by their spatial average, B and g, and assume 
the Abrikosov lattice for the spatial variation of the gap 
function, we obtain the analytic solution of the quasiclas- 
sical Green function. Using h 
the result is given by^^ 
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where i is the band index (cither a or tt), lj^ is the 
fermionic Matsubara frequency, W{z) — e^^ erfc(— iz) 
is the Faddeeva function, and Ai is the spatial average 
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of the energy gap in the £-band. We have assumed the 
isotropic s-wave gap for both bands. In eq.(l) we have 
introduced the quantity 

Vi±{k)^/2\e\B 

where k = kp/lkpl is the normahzed Fermi wave vector, 
and Vf±{k) is related to the Fermi velocity perpendicular 
to the applied field H. We shall give explicit expression 
of ve± (fe) later. 

The Maxwell equation and the gap equation follow 
from the stationary condition of the free energy with re- 
spect to B and A(, respectively. In the BPT approxima- 
tion the free energy measured from that of the normal 
state in the clean limit^^ is given by 



oe 
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where 



h = , , _ —AjWitUne), 



(3) 



(4) 



and (•••)£ denotes the average over the Fermi surface. 

In cq.(3), No£ is the DOS in the normal state, and 
7 ~ 0.5772 is the Euler's constant. The cut-off frequency 
Uc = QOttTc is used in the present paper. We have as- 
sumed B ^ H since MgB2 is the strongly type-II super- 
conductor with K ~ 26.^^ 

For simplicity, we use the cylindrical (spherical) Fermi 
surface for the two-dimensional (three-dimensional) a 
(tt) band. Then the Fermi velocities arc given by 

Va{k) = (vcTMb cos (f),v^^ab sin (f),v^^cSgn{kc)), (5) 

(fc) = f TT (sin 9 cos 0, sin Osincj), cos 9) , (6) 

where <p (9) is the azimuthal (polar) angle of k. Here the 
quasi two-dimensionality in a band is taken into account 
simply by Wcr,cSgn(A;c) with Vcr,c ^ v^^ab- We introduce 
ve±{k) by 



for H 



c, and 



(7) 



"v^+Xe^^iy, (8) 

for H J- c. Here xe characterizes the anisotropy of the 
coherence length perpendicular to the field. Thus, Xe = i- 
for II c. For H J. c there is strong anisotropy in cr 
band due to the quasi two-dimensionality. In this case, 
we determine xe = {£,ez / S.iy)'^ variationally by minimizing 
the free energy. ^^-^^ In the case of the Ginzburg-Landau 
theory, xeCL = (veJ/ivey)- 

Due to the coupling between a and tt bands, it may 
be expected that Xv differs from unity even if tt band is 
isotropic. It turns out, however, that the variational so- 
lution has no anisotropy in Xtt within the numerical ac- 
curacy, whereas Xo- has a strong anisotropy. This result is 



consistent with that obtained by Ichioka et ai, in which 
the local DOS(LDOS) of the tt band is nearly isotropic 
while that of the a band has strong anisotropy.^" On 
the experimental side, almost isotropic vortex core in 
the TT band is observed by scanning tunneling spec- 
troscopy(STS) at low fields, even in the case of i? _L c.^^ 
On the other hand, the use of the Abrikosov lattice with 
different xe's in two bands results in unphysical situ- 
ation since the structure of penetrating flux differs in 
two bands. However, after spatial average over a larger 
scale, this deficiency of the solution may not be serious 
as far as we are concerned with homogeneous properties. 
Since only spatially averaged quantities have a meaning 
in the BPT approximation, we use the variational solu- 
tion, = 1, throughout this paper. 

3. Anisotropy of thermodynamic quantities 

Let us first determine the strength of the pairing in- 
teraction in the matrix. 



V 



V 
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by fitting the specific heat data at = 0. Note that the 
present approximation recovers the ordinary two-band 
BCS theory^i'Si at i? 0. We use iVoa/A^o = 0.45 and 
No-jt/No = 0.55 (A^o is the total DOS) from the specific 
heat measurement.^ These values are in good agreement 
with the band structure calculation.^* Figure 1 shows 
the best fitting of the T dependence of the specific heat 
a.t H = 0. The coupling constants are determined to be 
N^V^ = 0.30, = 0.14 and V/V^ = 0.34. The solid 

circle is taken from the experimental data by Wang et 
al.^ We obtain the gap magnitude at T = as Aq-o = 6.8 
meV and A^ro = 2.4 meV, which arc in good agreement 
with those determined by tunneling spectroscopy.^'^ Our 
coupling constants are slightly different from those deter- 
mined by the band structure calculations'^^' ^'^ due to the 
simplification of the Fermi surfaces in the present study. 
In this work we try to reproduce of overall T dependence 
of the gap magnitude, i.e., specific heat, rather than de- 
termining the coupling constants uniquely. In fact, dif- 
ferent choice of coupling constants may give similar be- 
haviors in C(T).i5'29 

Next, we determine the Fermi velocities of two bands 
by fitting the anisotropy of the upper critical fields, 
r(r). The solid line in Fig. 2 represents our result for 
Va,clva,ab = 0.13 and VT^jvc^ab = 3.0, while the solid cir- 
cles are experimental data taken from Lyard et al}^ The 
result for Va-^c/va.ab = 0.13, UTr/Uff.ab = 1-0 is also shown 
in Fig. 2 (dashed line) for comparison. The slight modi- 
fication of the velocity ratio gives a considerable change 
in r(T) especially at high temperatures. The present 
result for T{T) is consistent with previous theoretical 
works. ^^"^^ It is emphasized that the strong tempera- 
ture dependence of r(T) originates from the interplay 
of the two gaps since the T dependence becomes much 
weaker in the absence of the interband coupling, V = 0. 
Note that r(r) is independent of T in the single-band 
anisotropic GL theory. The values of Va,c/'Va.ab — 0.13 
and v-j^jv^^ah = 3.0 are in good agreement with those 
obtained from the band structure calculation.^''' 
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Fig. 3. The in-plane field dependence of (a) the ZEDOS, and (b) the two gaps at T/Tc = 0.05. The solid circle represents the 
experimental data taken from Bouquet et al.^ For comparison the results for f jr/'i'CT.ai) = 10 are shown in (c) for ZEDOS and (d) for 
the gaps. 
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Fig. 1. The temperature dependence of the specific heat at H = 
0. The solid circle represents the experimental data of Wang et 
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Fig. 2. The ratio F of the upper critical fields Hc2(T) along the 
c-axis and the afe-plane against temperature. The solid circles are 
the experimental data by Lyard et al.^^ 



With these preliminaries, we discuss the field depen- 
dence using the coupling constants and the Fermi veloc- 
ities determined above. Figure 3 shows the in-plane field 
dependence of two gaps and the ZEDOS at T/Tc = 0.05. 
The solid circle denotes the experimental data taken from 
Bouquet et al.^ Owing to the smallness of A^r as com- 
pared with Act, the coherence length ^.n- in tt band is 
much larger than ^o-- Indeed, we obtain £,Tr/^a,ab ~ 7 and 
Ctt/Cct.c 60 for VTr/v„^ab — 3.0 usiug the BCS formula. 
On the other hand, the upper critical field is determined 
mostly by the pairing in the a band. Thus, the rapid rise 
of the ZEDOS should result as shown in Fig. 3(a). The 
contribution of the tt band to the ZEDOS saturates al- 
ready at H/H^^'^ ^ 0.05. Note that Ajr remains nonzero 
up to Hc2 because of the interband coupling as shown in 
Fig. 3(b). 

The results in the case of v-^jvc^ah — 1-0 are shown 
in Fig. 3(c) and (d) for comparison. The change of the 
velocity ratio only affects the low-field behaviors of the 
ZEDOS and of the gaps. On the contrary, r(T) of the iJc2 
ratio is affected much in higher temperatures as shown in 
Fig. 2. Experimentally Gonneli et al?^ report the rapid 



suppression of the gaps, while Bugolavsky et al. report 
no considerable suppression at low fields. There seems 
sample dependence in S,t! /S,a presumably due to the im- 
purities. There are a number of investigations on the field 
dependence of the ZEDOS and the gaps both theoreti- 
cally^^"^*' and experimentally.'^^' '^^ All of these show sim- 
ilar behaviors except in low-field region. 

Once we determine the necessary parameters, we cal- 
culate various quantities under magnetic fields without 
any free parameters. The T dependences of the specific 
heat under magnetic fields with v-^jva^ah = 3.0 are shown 
in Fig. 4 for (a) || c and (b) for H '_L c. For both field 
directions, the overall behavior of C{T,H) is consistent 
with the experimental results indicated in the figures. In 
particular, the observed specific heat jump at Tc{H) is 
reproduced very well for both field directions. 
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Fig. 4. The T dependence of the specific heat under magnetic 
fields (a) for H \\ c and (b) for H A- c. The experimental data 
represented by symbols are taken from Lyard et al}^ The exper- 



imental value of Hc2 at T = is H 
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Finally, we discuss the H and T dependence of Xa for 
in-plane fields. The T dependence of Xo- at = 11^2'^^ is 
shown in Fig. 5(a), where Xgl ~ 

V^V^,c/Va,ab ^ 0.18. At 

low temperatures the vortex cores elongate along c-axis 
slightly more as compared with the GL theory. In the 
limit of T Tc, the anisotropy recovers the GL value. 
The vortex cores also elongate monotonously along c-axis 
with increase of the field as shown in Fig. 5(b). It should 
be noted that the T and H dependence of Xo- does not 
come from the effect of two-band superconductivity, but 
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Fig. 5. (a) The temperature dependence of xir''^ at H = h'^^K 

1/2 

(b) The field dependence of x<t at T/Tc = 0.05. The anisotropy 

1 /2 

in the GL theory is given by Xql • 



comes purely from the quasi two-dimensionality of the a 
band. This conclusion follows since there is no essential 
difference even in the absence of the inter-band interac- 
tion, V = 0. 

4. Summary 

We have investigated the thermodynamic properties 
and its anisotropy of the two-band superconductivity in 
MgBj on the basis of the approximate analytic solution 
in the quasiclassical formalism. Using the strength of the 
pairing interactions and the Fermi velocities of two bands 
as fitting parameters, we have reproduced the T- and H- 
dependence of the ZEDOS, the two gaps, the specific 
heat, and the anisotropy of the upper critical field in 
wide regions of the phase diagram. 

The change of the velocity ratio affects considerably 
the low-field behaviors of the ZEDOS and the energy 
gaps. At the same time the velocity ratio determines the 
high-temperature behavior of the Hc2 anisotropy r(T). 
Thus, our theory predicts a strong correlation between 
the low-field behaviors of the ZEDOS and the high- 
temperature behavior of r(T). We have shown that the 
anisotropy of the coherence length in the a band has H 
and T dependence different from that in the GL theory. 
The isotropic nature of the coherence length in the tt 
band remains even if the inter-band coupling is present. 
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